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1. Introduction distinguish radio 'active' from radio inactive dwarfs. Possible 

physical characteristics include mass, temperature, activity, and 
(N In the past few years, a number of very low mass stars and brown ^^^^^^^^ ^^^^ Intriguingly, all active dwarfs have been found to 
T-H ; dwarfs (collectively termed ultracool dwarfs) have been con- ^^^^ ^^^^ ^ ■ ^^^^^^ j^^^ ^^gg^^j^ possibiUties, a depen- 
^ . firmed as a new class of radio active objects. A surpnsing feature ^^^^^ ^^^^^^^^ ^^j^^i^y^ ^ dependence on inclination angle. 
. ^ of these observations is the detection of penodic pulses of 100% ^he former suggests that slower rotators have weaker dynamo 
X ■ circulai-ly polarized emission (Hallmmetal, 2006, 2007, 200& ^^^j^n and hence weaker magnetic fields of insufficient strength 
H ■ l Berseretal. 1 |2009al | Doyleetal. | |2010t [ McLean et al. | |2011|). to produce detectable radio emission at the required frequency. 
. ^ .■ These penodic pulses have been confirmed to be produced by the ^^^^^^ ^^^^^^^ jj^pji^^ ^ dependence on inclination angle, 
electron cyclotron maser (ECM) instability, the same mechanism j ^ geometrical selection effect is associated with the highly 
that IS known to produce the planetaiy radio emission at kHz and ^^^^^^ ^^^^^ emission. A recent case study of three pulsing 
MHz frequencies (Treumami2QQ6), but it requires much more ultracool dwarfs does indeed confirm that all three have very 
powerful kilogauss magnetic fields. Unpolanzed and seemingly j^igh values of inclination (greater than 65 degrees). However, 
quiescent radio emission was present in the observations of all ^jj ^^^^^ ^^^f^ ^j^^ confi rmed to be very rapid rotators w ith 
detected dwarfs and has alternatively been attnbuted to gyrosyn- ^^^^^^^ ^^^^^^^^ < 3 ^^^^^ dHalfinan et al.ll200^ . 
chrotron emission and the ECM emission (see the above refer- 
ences). Clarifying the relationship between v sin ; and radio activ- 
Although much progress has been made in understanding the ity is imperative. If a dependence on rapid rotation underHes 
nature of the pulsed radio emission from these dwarfs by using the observed correlation between radio luminosity and v sin /, 
its diagnostic potential, it remains unclear which chai-acteristics this implies that the rotation-activity relationship, which is well- 

established for main-sequence stars, extends into the substellar 

Send offprint requests to: A. Antonova, tony@phys.uni-sofia.bg regime. If, on the other hand, the correlation between radio lu- 
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ABSTRACT 

Aims. We aim to increase the sample of ultracool dwarfs studied in the radio domain to allow a more statistically significant under- 
standing of the physical conditions associated with these magnetically active objects. 

Metliods. We conducted a volume-limited survey at 4.9 GHz of 32 nearby ultracool dwarfs with spectral types covering the range M7 
- T8. A statistical analysis was performed on the combined data from the present survey and previous radio observations of ultracool 
dwarfs. 

Results. Whilst no radio emission was detected from any of the targets, significant upper limits were placed on the radio luminosities 
that are below the luminosities of previously detected ultracool dwarfs. Combining our results with those from the literature gives 
a detection rate for dwarfs in the spectral range M7 - L3.5 of ~ 9%. In comparison, only one dwarf later than L3.5 is detected in 
53 observations. We report the observed detection rate as a function of spectral type, and the number distribution of the dwarfs as a 
function of spectral type and rotation velocity. 

Conclusions. The radio observations to date point to a drop in the detection rate toward the ultracool dwarfs. However, the emission 
levels of detected ultracool dwarfs are comparable to those of earlier type active M dwarfs, which may imply that a mildly relativistic 
electron beam or a strong magnetic field can exist in ultracool dwarfs. Fast rotation may be a sufficient condition to produce magnetic 
fields strengths of several hundreds Gauss to several kilo Gauss, as suggested by the data for the active ultracool dwarfs with known 
rotation rates. A possible reason for the non-detection of radio emission from some dwarfs is that maybe the centrifugal acceleration 
mechanism in these dwarfs is weak (due to a low rotation rate) and thus cannot provide the necessary density and/or energy of ac- 
celerated electrons. An alternative explanation could be long-term variability, as is the case for several ultracool dwarfs whose radio 
emission varies considerably over long periods with emission levels dropping below the detection limit in some instances. 

Key words. Stars: low-mass, brown dwarf - Radio continuum: stars ~ Radiation mechanism: general ~ Stars: activity 



2 



A. Antonova et al.: A Volume Limited Radio Survey of Ultracool Dwarfs 



minosity and v sin / indicates a geometrical selection effect, this 
implies that very strong magnetic fields (kG) are ubiquitous in 
the substellar regime, independent of rotation rate. Thus, to clar- 
ify the relationship between v sin / and radio activity, a larger 
more statistically significant sample of pulsing dwarfs must be 
established. 

The ECM emission has proved a vital diagnostic tool for re- 
mote sensing of the magnetic field strength and topologies of 
ultracool dwarfs. The ECM emission is generated at the elec- 
tron cyclotron frequency denoted by Vc » 2.8 x 10^ B Hz (if the 
electrons are not relativistic), enabling measurement of the mag- 
netic field strength in the source region of the pulsed emission. 
This led to the realiza tion of kG magnetic fields in late-M dwarfs 
(iHallinan et alJl2007l) and subsequently the first confirmation of 
kG magnetic fields for an L dwarf, the latter establishing strong 
magnetic dynamo action out to spectral type L3.5 ( Hallin an et al.l 
l2008h . 

Very recently, iRoute & WolszczanI (1201 2l) reported the first 
detection of radio emission from a dwarf of spectral type later 
than L3.5 - the T6.5 brown dwaif 2MASS J1047539+212423. 
They have detected circularly polarized bursts at 4.75 GHz with 
the Arecibo telescope and invoked the ECM mechanism as the 
most likely source of the emission. This latest observation con- 
firms that detecting ECM emission remains the most promising 
method to measure magnetic field strengths in cooler late-type 
L and T dwarfs. Reiners & Basri ( 2007.) measured the magnetic 
field strengths of a number of late-M dwarfs through the mea- 
surement of the Zeeman broadening of the magnetically sensi- 
tive Wing-Ford FeH band. However, this technique encounters 
difficulties when applied to L and T dwarfs due to the heavy 
saturation of the FeH lines. Therefore, the continued search for 
radio emission from late-L and T dwarfs is essential to the efforts 
to diagnose the strength of magnetic fields in these objects. 

To address the above questions we conducted a volume- 
limited survey of 32 ultracool dwarfs of spectral types M7 to 
T8. 



2. Observations and data reduction 

The target list was compil ed using the 2MU2 20 pc volume lim- 
ited catalog of ICruz et al. I (|2007). We selected field dwarfs that 
had no previous observations at radio frequencies and whose 
declination is > -30° (the latter ensuring that the observing time 
is sufficiently long for achieving the goals described below). This 
yielded 14 dwarfs in the spectral range M7 - L3.5, with a view 
to expanding the current sample of detected pulsating ultracool 
dwarfs. We also selected 18 dwarfs of types L4 to T8 with the 
aim to detect radio emission from a source later than L4. Their 
properties are listed in Table [1] 

The observations were conducted with the NRAO Very 
Large ArrajQ in the period 09 August - 09 September 2009 us- 
ing the standard continuum mode with 2 x 50 MHz contigu- 
ous bands. The observing frequency was chosen to be 4.9 GHz, 
because any ultracool dwarf observed at both 4.9 GHz and 8.5 
GHz has been found to produce quiescent emission at 4.9 GHz 
of the same luminosity or higher than that detected at 8.5 GHz. 
Furthermore, in the case of the ECM emission, an upper cut-off 
frequency is expected for the radio emission that is dependent 
on the maximum magnetic field of the dwarf. A survey at 4.9 
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GHz thus allows for a lower magnetic field strength, which ad- 
ditionally increases the chance of detection. For each source, a 
suitable phase calibrator was chosen from the VLA calibrator 
manual while the flux density scale was determined using the 
calibrators 3C48, 3C138, and 3C286. 

We observed each target in the survey for a duration deter- 
mined by its distance, thereby ensuring a common upper limit 
on the luminosity of the radio emission from each dwarf. For the 
dwarfs spanning the spectral range M7 - L3.5, the duration was 
chosen to match the one required to place an upper limit on the 
radio luminosity of ~ 10'^ ergs s"^ Hz"^ i.e., lower than each 
of the previously detected pulsing ultracool dwarfs. Considering 
that no detection has previously been reported for late L and T 
dwarfs, we chose to observe each of the later type (> L4) dwarfs 
for the duration required to place an upper limit on radio lumi- 
nosity of ~ 5 X 10'^ ergs s"' Hz"', i.e., to partially allow for any 
drop off in emission that might occur with later spectral type. 
All targets were observed in a dynamic scheduling mode with a 
time resolution of 3.3 sec. This mode was therefore sufficient to 
detect short radio bursts, in addition to any quiescent emission. 

Data reduction was carried out with the Astronomical Image 
Processing System (AIPS) software package using standard rou- 
tines. The visibility data were inspected for quality both before 
and after the calibration procedures, and noisy points were re- 
moved. To make the maps we used the task IMAGR and then 
CLEANed the region around each source in the individual fields. 

3. Results 

We found no sources at (or near) the positions of the stud- 
ied dwarfs in the respective maps, i.e., we did not detect ra- 
dio emission from any of the observed dwarfs. As discussed in 
Section |2] in the case of non-detections we aimed to obtain data 
with a sensitivity that would allow us to place upper limits of 
Ly < 10'^ ergs s ' Hz ' for the dwarfs in the range M7 - L3.5 
and Ly < 5 X 10'^ ergs s ' Hz ' for the later type dwarfs (L4 - 
T8). We find that for all dwarfs in the sample whose distances 
are within ~ 12 pc, the detection limits are at or below these val- 
ues. For the most distant dwarfs, however, the limits were not 
reached due to insufficient sensitivity. The 3cr upper limits on 
the flux and the respective radio luminosity limits are shown in 
TablelU 

We plot our results together with those from previous studies 
of UCDs at 4.9 GHz in Figure[T] This work increases the number 
of studied objects by more than 150% and lowers the upper lim- 
its on the radio luminosity. There is only one upper limit whose 
value is above the value of detected emission for a given spectral 
type - that for the M8.5 dwarf LSR J1826+3014, which is one of 
the most distant dwarfs i n the survey (d = 13.9 pc). 

iMcLean et al.l (l2012h recently reported on a 8.5 GHz radio 
survey of 100 M and L dwarfs spanning the range M4 to L4 (76 
of which are ultracool dwarfs) and the detection of three new 
radio-emitting late-type dwarfs. Comparing the target list of their 
observation to ours, we found seven dwarfs in common, none of 
which was detected at the higher frequency either. We list these 
dwarfs and their flux and luminosity upper limits in the respec- 
tive frequencies in Table|2] We have calculated the upper limit on 
the 8.5 GHz radio luminosit y using the distances l isted in Table 
12] and the fluxes reported bv iMcLean et al.l (l2012l) . Despite the 
small differences in the 3cr limits on the radio flux between the 
two surveys, we find that the limits in radio luminosities agree 
well. 

In addition. lMcLean et alj (1201 2|) compiled a list of all previ- 
ous radio observations of late-type dwarfs. Combining their list 
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Table 1. Sample of 32 dwarfs, their properties, and results from the present survey. 



Name 



Other name 



Sp. T. 



d v sin / 

(pc) (km 



Lboi 
(Lo) 



(Lh,t / Lbol) P (4.9 GHz) 



Lv,4.9 

(erg Hz-') 



< 1.53 XIO" 

< 1.51 xIO'^ 
<7.43 xIO'^ 

< 1.06xI0'2 
<9.73 xlO'2 

< 1.41 xlO'^ 

< 1.37 xIO'^ 

< 7.54x10'^ 
<2.0I xIO'3 

< 1.04x10'^ 

< 9.50x10'^ 

< 1.14x10" 

< 1.01 xlO'^ 

< 8.01 xlO'- 

< 1.04x10" 
<4.83 xlO'- 

< 8.55 xIO'- 

< 9.86x10'- 
<4.82xl0'2 
<4.70xI0'2 

< 3.36 xlO'^ 

< 7.38x10'^ 

< 1.03 xlO'^ 

< 1.78 xm" 

< 3.68 xIO'2 

< 6.14x10'^ 

< 5.28x10'^ 

< 2.98x10'^ 
<2.54xI0'3 
<5.34xl0'2 

< 3.52x10'^ 

< 7.18x10'^ 

< 4.89x10'^ 

< 1.84x10'^ 
<2.23 xlO'^ 



2MASS J 1 048 1 258- 1 1 20082 
2MASS J 1757 1539+7042011 
2MASS J11554286-2224586 
2MASS J05395200-0059019 
2MASS J12505265-2121136 
2MASS J04351455-1414468 
2MASS J02 150802-30400 11 
2MASS J05392474+4038437 
2MASS J18261131+3014201 
2MASS J14284323+3310391 
2MASS J17312974+2721233 
2MASS J0921 1410-2104446 
2MASS J08283419-1309198 
2MASSI J0700366+3 15726 
2MASS J05002100+0330501 
2MASS J04351455-1414468 
2MASS J()3552337+l 133437 
2MASS J05395200-0059019 
2MASS J17502484-0016151 
2MASS J21481633+4003594 
2MASS J09083803+5032088 
2MASS J08575849+5708514 
2MASS J02572581-3105523 
2MASS J08300825+4828482 
2MASS JO 1365662+0933473 

2MASS J15031961+2525196 
2MASS J09373487+2931409 
2MASS J21543318+5942187 
2MASS J07271824+1710012 
2MASSJ1 1145133-2618235 

2MASS J09393548-2448279 

2MASS J14571496-2121477 



GJ 3622 
LP 44-162 
LP 851-346 

SDSS J05395 1.99-005902.0 
DENIS-P J125052.6-212113 

LHS 1367, LP885-35 
LSR J0539+4038 
LSR J 1826+30 14 
GJ 3849, LHS 2924 
LSPM J1731+2721 
DENIS-P J0921 14. 1-210445 
DENIS-P J082834.3-130919 



SDSS J05395 1.99-005902.0 



SDSS J085758.45+570851.4 
DENIS-P J025725.7-3 10552 
SDSS J083008. 12+482847.4 
SIMP J013656.5+093347.3 



GL570D 



M7 

M7.5 

M7.5 

M7.5 

M7.5 

M8 

M8 

M8 

M8.5 

M9 

LO 

L2 

L2 

L3.5+L6 

L4 

L4.5 

L5 

L5 

L5.5 

L6.5 

L7 

L8 

L8 

L9 

T2.5 

T5.5 

T6 

T6 

T7 

T7.5 

T8 

T8 



4.5 

12.5 

9.7 

3.84 

11. 1 

14 

12.37 

10 

13.9 

11.8 

11.8 

12 

11.6 

12.2 

13.03 

9.8 

12.6 

13.1 

8 

7.9 

15.9 

11 

9.6 

13.0 

6.4 

8.27 

7 

6.14 

18.8 

9.09 

7 

10 

8.7 

5.34 

5.88 



3.0 

33 
33 



10 
15 
15 
33 
41 



10 

34 



31 



36 



-4.63 
-5.01 
-4.58 



-4.7 
-4.6 
< - 6.42 
-5.68 



30 



-3.16 
-3.48 



-3.55 



-3.62 
-3.74 
-4.01 

-3.96 
-4.26 
-4.12 
-4.03 
-4.2 

-4.07 



-4.82 
-4.58 



- 5.28 

-5.26 <-5.4 

- 5.69 

-5.52 - 



<6.3 
<8.1 
<6.6 
<6.0 
<6.6 

< 6.0 
<7.5 
<6.3 
<8.7 
<6.3 
<5.7 
<7.2 
<6.3 
<4.2 
<5.1 
<4.2 
<4.5 
<4.8 
<4.3 
<6.3 

< II. 1 
<5.1 
<9.3 
<8.7 

< 7.5 

< 9.0 
<6.6 
<6.0 
<5.4 

<6 

<5.4 
<5.4 



The columns are (left to right): name of the object; other name; spectral type; distance; projected rotational velocity; bolo- 
metric luminosity; Ho- luminosity; radio flux at 4.9 GHz; radio luminosity at 4.9 GHz. Sp ectral types, distances and v sin i 
measurements taken from: SIMBAD, Bemat et al. (2010); Goldman et al. (|2010|);|Lee et al. (2010); Reiners & Basri' ( l2010h; 
ISeifahrt et al. (2010); Smart et al. (2010); Del Burgo et al. (2009); Faherty et al. (2009); Leggett et al. (2009); Stephens et^ 
([2009); Looper et al. (2008); Reid et al. (20 08);,Reine rs & Basri ( 2008); Blake et al. ( 2 007) ; Kendall et al . ( 2007 ) ; Looper et^ 




)8);, Reine rs 

Reiners & B asri (2007); Schmidt et al.l (1200 7); Artig au et al. 

Zap atero Osorio et al 1 12006); iCrifo et al 



R eid et al 
iTinney et all ( 120051 
120031120021) 



112006), 

I); iGolimow ski et al.l ( l2004h : IVrba et al.l ( 12004ft ; iBurgasser et al.l ( l2003al) ; ICmz'et a l 




2006); 
2005.); 



Hen ry et al. 1 2006h; iMora les-Calderon et al 



.Ctishing et alJ (1200: 



'Lodieu et aL 42005.. 
(2003); lL6'pine etal.l 



with the present survey raises the total number of studied ultra- 
cool dwarfs to 193, with 13 detected to posses radio emission 
(see Figure |4|i. Although this might suggest a detection rate of 
~ 6% in the ultracool dwarf regime, a more detailed analysis in- 
dicates otherwise. All but one of the detected ultracool dwarfs 
span the spectral type range M7 - L3.5. Therefore, the detection 
rate over this spectral range is 12 confirmed sources from 134 
observations, i.e., ~ 9%. In comparison, only one dwarf later 
than L3.5 was detected in 53 observations. 

In Fig.|2] we plot the number distribution of cool dwarfs later 
than MO (red line) observed in the radio domain compared to 
those with a detected radio signal (blue line) with respect to their 
spectral type. We include earlier type M dwarfs to show the over- 
all trend in detection rate. The insert panel in the figure shows 
the detection rate as a function of spectral type, defined by Nr/N, 
where A^i is the number of dwarfs with detected radio signal in a 
certain spectral type range; is the number of observed dwarfs 
in the same spectral type range. To estimate the uncertainty of 
the detection rate as a function of spectral type, we adopted the 
binomial distribution with the limits equivalent to Icr hmits for 
a Gauss ian distribution (see, e.g., the article of IBurgasser et al.l 
(l2003bl) . where this model was applied to the statistics of T- 
dwarf binaries). There is a drop in detection rate from ~ 50% 



at spectral type M4 to less than 15% for later spectral types. 
The drop may be due to the small sample of observed UCDs, 
especially those later than LO. Alternatively, the drop may have 
physical implications. We discuss these possible implications in 

The top panel in Fig.[3] shows the observed radio luminosity 
as a function of rotation velocity for the M and L dwarfs (most 
of them are observed by the VLA, mainly at 4.9 and 8.5 GHz). 
Solid circles are for the M and L dwarfs with detected strong 
radio signal, whilst open triangles represent the upper-limit of 
the luminosity of some M and L dwarfs. The relation between 
the luminosity and v sin / in this figure is not very clear. It seems 
that the number of M dwarfs possessing strong radio emission 
is more than the number of L dwarfs. Even for one individual 
dwarf, the detected radio luminosity can vary over a wide range 
(about 1-2 orders of magnitude) in a short observation period, 
which may suggest that the radio emission is strongly associated 
with its generation mechanism and may not be directly corre- 
lated with rotation rate. 

The bottom panel in Fig. [3] illustrates the number distribu- 
tion of all observed M (red line) and L (blue line) dwarfs as a 
function of rotation velocity (v sin /). The number distribution of 
M dwarfs decreases with respect to v sin /, perhaps caused by 
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Fig. 1. Radio luminosity - spectral type plot for all ultracool dwarfs observed at 4.9 GHz. Filled triangles show the upper limits on 
radio luminosities of the dwarfs from the present survey, open tri angles mark upper limits from the literature, and filled squares 
represent radio luminosities of detected dwarfs from the literature (lBergerll2006i iBurgasser & Putmanll2005[ lAntonova et al.ll2007l 
[2008; Route & Wolszczan 2012, and references therein). 

Table 2. Upper limits on the radio flux and radio luminosity at 4.9 GHz (present survey) and 8.5 GHz (iMcLean et al.ll2012l) of the 
seven dwarfs that are present in both surveys. We used the distances listed in Tabl^Uto calculate Ly g j. 



Name 


Other name 


F(4.9 GHz) 

(10-^Jy) 


Ly,4.t) 

(erg s"' Hz"') 


^(8.5 GHz) 

(10-^Jy) 


Lv,8.5 

(erg 


Hz-') 


2MASS J10481258-1120082 


GJ 


3622 


< 


6.3 


< 


1.53 xlO'^ 


< 


9.6 


< 


2.33 


xlO'^ 


2MASS J17571539+7042011 


LP 


44-162 


< 


8.1 


< 


1.51 xlO" 


< 


11.7 


< 


2.19 


xlO'3 


2MASS J11554286-2224586 


LP 


851-346 


< 


6.6 


< 


7.43 xlO'2 


< 


9.0 


< 


1.01 


xlO'3 


2MASS J12505265-2121136 






< 


6.6 


< 


9.73 xlO'2 


< 


7.2 


< 


1.06 


xlO'^ 


2MASS J09211410-2104446 






< 


7.2 


< 


1.14x10'^ 


< 


7.5 


< 


1.29 


xlO'^ 


2MASS J08283419-1309198 






< 


6.3 


< 


1.01 xlO'3 


< 


6.6 


< 


1.06 


xlO'3 


2MASSI J0700366+3 15726 






< 


4.2 


< 


8.01 xlO'^ 


< 


7.8 


< 


1.39 


xlO'3 



a stellar-wind-assisted magnetic braking. More observations are 
needed to exclude this selection effect. Also, since the major- 
ity of the data represent upper limits, no statistically significant 
conclusions can be inferred at present. 



4. Discussion 

In this section, we discuss the physical implication of our survey, 
including the influence of the radiation mechanism, the geomet- 
ric effect, and the rotation of the dwarfs. A brief discussion on 



two other activity characteristics - X-ray and Ha emission from 
the radio observed sample - is also included in Appendix lAl 

4.1. Radiation meclianisms and geometric effect 

For 12 of the 13 radio active dwarfs a quiescent component of 
the radio emission was also detected. The T6.5 dwarf 2MASS 
J1047539-I-212423 was observed with the Arecibo single-dish, 
which has a low sensitivity to unpolaiized quiescent emission, 
thus the presence of such a component is not confirmed, but it 
is likely. The nature of the quiescent emission is still debated. 
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Fig. 2. Number distribution of observed dwarfs (red line) and dwarfs with detected radio signal (blue line) with respect to their 
spectral type. The insert panel shows the detection rate as a function of spectral type, defined by Nr/N, where is the number of 
dwarfs with detected radio signal in a certain spectral type range; is the number of observed dwarfs in the same spectral type 
range. In the spectral type range of M0-T8, the total number of observed dwarfs is 273, while the number of dwarfs with detected 
radio signal is 34. The zero detection rate around spectral type M2 is most likely due to the limited number of observations. Also 
shown are the uncertainties of the detection rates, estimated by adopting the binomial distribution at Icr level for the Gaussian 
distribution (see §[3]i. 



It could either be due to an incoherent process such, as gy- 
rosynchrotron radiation, or depolarized electron cyclotron maser 
instability emission (Hallinan et al. 2008). On the other hand, 
the ECM is confirmed to be the dominant mechanism driv- 
ing the bulk of the emissio n from five of the detected dwarfs 
dHallinan et alj |2008: BerHeretalj|2009b; McLean et al. 2011). 
For the L3.5 dwarf 2MASS J0036 1 6 1 7+ 1 82 1104, H alhnan et al J 
(|2QH) showed that the radio emission is originally due to the 
ECM, but perhaps depolarized during propagation toward the 
observe r. In addi t ion, re cent particle-in-cell numerical simula- 
tions by Yu et al.l (1201 2h have shown that, under certain circum- 
stances, maser emission can be intrinsically weakly polarized. 
The same mechanism can be applied to generate the bright, 
100% polarized burst/flare detect ed from another ultracool d warf 
- DENIS-P J104814.9-395604(Burg asser& Putmanl l200l. 

Particle-in-cell simulations, which are based on the particle- 
field theory, indicate that a mildly relativistic or relativistic 
electron-beam-driven cyclotron maser can be an effective mech- 



anism to release electrom agnetic waves and heat the surrounding 
plasmas (lYu et al.ll2012l') . Furthermore, from the diffusion pro- 
cess of the electrons in velocity space, a high-energy tail may 
be rapidly developed along the direction perpendicular to the 
magnetic field, which can eventually evolve to moderately or 
strongly relativistic electrons depending on the initial energy of 
the electron current, and contribute gyrosynchrotron radiation. 
This may lead to the appearance of a radio continuum and the 
deformation of the spectral energy distribution. Also, the sim- 
ulations show that a series of discrete spectral lines can appear 
at certain frequency bands, which may increase the difficulty of 
finding the fundamental cyclotron frequency in the observed ra- 
dio frequencies. 

However, it is likely that the magnetic field inhomogeneity 
will smooth out these bands, thus producing a continuous spec- 
trum. Therefore, it is possible that determining the cyclotron fre- 
quency and detecting radio emission from a larger number of 
UCDs may require observations over a wider frequency band. 
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Fig. 3. Top panel: Radio luminosity as a function of rotation velocity for the UCDs observed in the radio domain (mainly at 4.9 
and 8.5 GHz). Red solid circles are for M dwarfs, the blue solid circles represent L dwarfs. The solid circles connected by black 
lines denote different measurements for the same target. Red open triangles denote the upper-limit of the luminosity o f M dwarfs; 
blue o pen triangles stand for the upper-limit of the luminosity of L dwarfs. The data are mainly taken from Table[T] McLean et al.l 
(l2012h . and references therein. Bottom panel: the number distribution of the observed dwarfs as a function of rotation velocity. The 
red solid line is for M dwarfs (normalized by the total number 91), the blue dashed line is for L dwarfs (normalized by the total 
number 48). 



In the case of beamed emission, the geometric selection may 
play an important role in determining the detection rate (see in- 
sert panel in Fig. |2]i. The probability P to detect this emission 
can be expressed as P = 26 /(4n). The half solid angle of the 
radiation cone 6 can be approximated by cos 6 » j, where 



and c are the velocity of electrons in a local region and the speed 
of light. The current detection rate for dwarfs over spectral type 
M7 - L3.5 is about ~ 9%, which may imply that we need a 
mildly relativistic electron beam with an average velocity of at 
least 0.8c. 
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Fig. 4. Radio luminosity - spectral type plot for all ultracool dwarfs observed at radio frequencies. Filled triangles show the upper 
limits on radio luminosities of the dwarfs from the present survey. Open t riangles mark upper limits from the literature and filled 
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The numerical simulations indicate that a relativistic elec- 
tron beam can release more than 30% of its kinetic energy to 
electromagnetic energy in a rela tively weak magnetic field (hun- 
dreds of Gauss) ( lYu etal.1120121) . So one tentative conclusion is 
that fast electron beams may be common even in cool objects, 
e.g., L dwarfs. The electrons may come from the magnetic field 
coupled ionizat ion (or discharging) process in the atmosphere, 
as discussed in iHelling et al.l (1201 ih . Alternatively, they can be 
associated with the internal activity of the ultracool dwarfs. 

The ECM emission and/or geometric selection may be the 
explanation for the lack of detection of two particular dwarfs 
in our list - the M9 dwarf LHS 2924 and the M7 dwarf GJ 
3622, both listed in SIMBAD as flare stars. LHS 2924 has de- 
tected Ha emission (see Table [U, magnetic field strength of Bf 
= 1.6 + 0.2 kG (where f is the fi lling factor), and v sin / - 
10 km s ' (iReiners & Basrill2007h . Depending on the inclina- 
tion angle, the rotational velocity could be in the range 10 to 57 
km s ' . Considering the presence of both chromospheric activity 
and kilogauss fields, it would not be surprising if radio emission 
was present as well. However, this is not the case. We have not 
detected this dwarf down to a 3<t flux level of 6.3 xlO"^ Jy at 
4.9 GHz. lMcLean et alJ ( l2012l) reported a 3cr upper limit of 8.4 
xlO"^ Jy at 8.5 GHz. The lack of detection at the higher fre- 
quency is not surprising considering the magnetic field strength 
measured. Detection at 4.9 GHz, on the other hand, should have 



been possible since the surface magnetic field is not less than 
1600 G (i.e., the cyclotron frequency exceeds 4.5 GHz). Thus, 
the absence of maser emission may be attributed to geometric ef- 
fects since the beamed radio emission might never sweep in the 
direction of Earth during a full rotation of the UCD. However, 
for TVLM 513 the pulses have a duty cycle of less than 15%, 
thus as an alternative explanation we may have missed the range 
of rotational phase during which pulsed emission is present. 

The second dwarf, GJ 3622, has a detected Ha emission (see 
Table [TJ and X-ray emission (log(Lx) - 25.96 erg s"'), a mag- 
netic fieldstrength of Bf - 600 + 200 G, and v sin / = 3 km s ' 
(Re iners & Ba sri 2010; Lee et al.l2010l:ISchmitt & Li efke 2004j). 
The absence of radio emission from GJ 3622 can be explained 
by the fact that the surface magnetic field is only about 600 G, 
so that the maximal cyclotron frequency is 1.7 GHz, which is 
below our observation frequency. 

The above considerations, however, cannot explain the ab- 
sence of gyrosynchrotron emission, which has a wide directiv- 
ity pattern and can be produced at frequencies well above the 
cyclotron frequency. Therefore, despite a strong magnetic field, 
the concentration and/or energy of accelerated electrons in the 
magnetospheres of the mentioned dwarfs may be insufficient to 
produce observable radio emission. 
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4.2. Rotation 

All radio-detected dwarfs have been found to have high v sin / 
(>20 km s ') valu es which may be int erpreted as a dependence 
on rapid rotation (McLean et al.ll2012l) . However, over 50 ultra- 
cool dwarfs with measured v sin i > 15 km s ' have not been 
detected as radio sources. It seems that fast rotation is a suffi- 
cient, but not necessary condition for late-M and L dwarfs to be 
radio-active. This tentative conclusion is consistent with the ro- 
tation and magnetic activity relation for a sample of M-dwarfs 
determined by measuring the equ ivalent widths of the Call H 
and K lines ( Browning et alJl2010l) . 

The radio emission from ultracool dwarfs can be strongly as- 
sociated with their rotation rate in two ways. First, rotation plays 
a key role in the stellar dynamo, thus determining the magnetic 
field strength. Second, rotation can directly affect the process of 
particle acceleration. For the planets of the solar system, the en- 
ergy of the accelerated particles can come from three sources: a) 
interaction of a magnetosphere with the solar wind (this process 
occurs at all planets and is especially important for the magne- 
tospheres of the Earth, Uranus, and Neptune); b) interaction of a 
magnetosphere with a satellite (like in the lo-Jupiter system); c) 
centrifugal acceleration of the magnetospheric plasma. 

The latter effect is caused by the fact that in a rigidly ro- 
tating magnetosphere, the centrifugal force exceeds gravity at 
some distance from the star/planet. This results in an accumula- 
tion of plasma and in an increase of the plasma pressure in those 
regions. Eventually, the plasma pressure exceeds the magnetic 
field pressure, so that the magnetic field lines become broken. 
As a result, the magnetic field at distances exceeding a certain 
radius does not rotate with the same rate as the star/planet; cur- 
rent sheets are then formed at the boundary between the inner 
(co-rotating) and outer parts of the magnetosphere. The electric 
fields in the mentioned current sheets, in turn, can accelerate 
electrons up to relativistic ener gies; thus the rotat i on energy is 
converted into particle energy (Andre et al. 19881 iLinskv et al.l 
[1992; Usov & Melrose 1992; Leto et al. 2006). 

The above process has been proven to be important for the 
magnetospheres of Jupiter and (to a lesser degree) Saturn. It is 
also proposed as the main particle acceleration mechanism for 
the radio-emitting star CU Virginis and other similar objects. 
Regarding the ultracool dwarfs, if a dwarf is not part of a bi- 
nary system, the interaction with an external stellar wind is obvi- 
ously absent. Satellite-induced radio emission can be recognized 
by its specific light curv es reflecting both the d warf's rotation 
and the satellite motion (iKuznetsov et aljr2012l) . The satellite- 
induced particle acceleration can be ruled out at least for some 
radio-emitting dwarfs, because this model cannot account for 
the strictly periodic emission pulses (which are typical of some 
dwarfs) nor for nearly constant quiescent emission (which is also 
observed). Thus, for an ultracool dwarf without close-in satel- 
lites, one probable particle acceleration mechanism is the cen- 
trifugal acceleration. 

Recently, iNichols et alJ (l2012h have estimated the parame- 
ters of the electric currents produced in the rapidly rotating mag- 
netospheres of ultracool dwarfs through angular velocity shear; 
it has been demonstrated that the power carried by these cur- 
rents can be sufficient to provide the observed intensities of radio 
emission. An assumption in this model is that the spin and mag- 
netic axes are co-aligned, which if the polarity were reversed, 
implies that only ~50% of UCDs would have auroral and radio 
emission. Coupling this wit h the known radio-a ctive duty cycle 
of ~ 0.15 phase for UCDs dHaUinan et al.ll2008l) . similar to the 



active duty cycle of Jupiter of 0.14 phase (iHiggins et al.lll996l) . 
could easily imply a very low detection rate. 

Evidently, the efficiency of this mechanism is strongly de- 
pendent on the rotation rate: an increase of the rot ation rate 
both increases the power transferred to the particles (ISchriiveii 
2009) and shifts the acceleration region toward the stellar sur- 
face (for the rotation period of a few hours, the co-rotation ra- 
dius is expected to equal just a few stellar radii). We propose the 
following explanation of non-detection of radio emission from 
the dwarfs LHS 2924 and GJ 3622: although the rotation rate of 
these dwarfs is relatively low, it is sufficient to produce the mag- 
netic fields of kilogauss strength (indeed, the dynamo models 
predict that the dynamo effects saturates at relatively low rota- 
tion rates corresponding to v about 5 km s so that a further 
increase in the rotation rate does not affect the magnetic field 
strength significantly (McLean et al. 2012)). At the same time, 
owing to the low rotation rate, the centrifugal acceleration mech- 
anism at these dwarfs is too weak and cannot provide the neces- 
sary density and/or energy of accelerated electrons. 

There is still no conclusive evidence on whether rotation or 
inclination angle is more important for detecting radio activity 
from UCDs. It might be the case that both play a significant role 
- the first for generating the emission and the latter for detect- 
ing it. The ultracool dwarfs with detected ECM emission are all 
found to be fast rotators with high inclination angles. The ques- 
tion whether detection of the radio emission from a larger sample 
depends on observing frequency may only be addressed via ob- 
servations of individual targets over wide frequency bands with 
the upgraded JVLA and ATCA. 

There are currently several ongoing modeling projects that 
aim at a better understanding of the emission process that pro- 
duces the periodic radio bursts with high brightness temperature 
and polarization degree in UCDs. Recent computations indicate 
that the rotation-modulated ECM emission can interpr et the ra- 
dio lig ht curve from die M8.5 dwarf TVLM 513-46546 (lYu et al.l 
1201 l l) therefore the properties of the radio emission region can 
be determined. Future work will focus on the combination be- 
tween different magnetic field topology and a particle simulation 
to reproduce the local radio emission region. 

IKuznetsov et alJ (1201 2l) presented results for the same dwarf 
where the emission properties were similar to the auroral radio 
emission of the magnetized planets of the solar system via the 
electron-cyclotron maser instability. Two models were consid- 
ered where the emission was caused by interaction with a satel- 
lite or derived from a narrow sector of active longitudes. It was 
found that the model of emission from an active sector is able 
to qualitatively reproduce the main features of the radio light 
curves; the magnetic dipole needs to be highly tilted (by about 
60°) with respect to the rotation axis. For the most often ob- 
servationally studied object, TVLM 513-46546, the model of 
the satellite-induced emission was inconsistent with the obser- 
vations. 

4.3. Other possibilities 

Another possible explanation for the lack of detected radio emis- 
sion could be long-term variability. To date, there are several ul- 
tracool dwarfs whose radio emission varies considerably over 
long periods, with e mission levels dropping below the detection 
limit in some cases dAntonova et al.l 120071 ; iBerger et al.l 120101; 
McLean et al. 2012i). 

Long-term variability, however, requires systematic monitor- 
ing of large sample of targets. Furthermore, full rotational phase 
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coverage should be obtained to account for the narrow beaming 
of electron-cyclotron maser emission. 
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Appendix A: Activity in j-lo' and X-rays of the sample of radio-observed dwarfs. 

In Table IATI we list all ultracool dwarfs observed in the radio domain, their distances, vsin i, radio fluxes, and luminosities/upper 
limits in Ha, X-rays, and the radio. About 60% of the radio-observed dwarfs also have observations in Ha, while only ^ 18% 
were observed in X-rays. Figures lATI and lAn show the Ha and X-ray quiescent luminosities as a function of spectral type. Despite 
the trend of lower hnal^hoi with later spectral type, only three UCDs of type < L7.5 have upper limits on their Ha luminosities. 
Toward the later spectral types, only one dwarf has a marginal Ha detection - the one detected with Arecibo, the T6.5 dwarf 2MASS 
J1047H-21 jRoute & Wolszczanll2012t iBurgasser et al.ll2003ah . However, its activity level is comparable to that of early-L dwarfs. 
Overall, these results point to chromospheric activity present throughout the L spectral type and possibly as far as late-T types. 

For late-type, main-sequence active stars (and especially active M dwarfs), the re is an empirical correlation between the 
observed X-ray and radio emission (Lr/Lx ~ 10 '^^) found bv lGuedel & Ben3 (Il993h . The suggestion is that the same population 
of relativistic electrons is responsible for both emissions - the spiraling motions of the electrons in the magnetic field give rise to 
incoherent radio emission prior to heating the coronal plasma, which leads to th ermal X-ray emissi on. However, for the ultracool 
dwarfs this correlation is found to be violated by several ord ers of m agnitude (Berger e t al.ll201(jl and references there in). For 
UCDs with observations in both the X-ray and radio domains, ISt elzer et al. (2012) suggested that they may be separated into two 
groups: X-ray flaring but radio-faint objects that are also slow rotators, and X-ray-faint but strong radio-bursting dwarfs with fast 
rotation. However, there are several dwarfs with vsin / in the range 15 - 60 km s ' that do not follo w this pattern. F or example, the 
L2 dwarf Kelu-1 (with vsin / - 60 km s ') has been reported to only possess X-ray emission (Au dard et alj|2007l) . As seen from 
Table lA.ll and Figure IaTSI the sample of dwarfs observed both in radio and X-rays is too small to allow any significant conclusions. 



Table A.l. Activity characteristics of all ultracool dvyarf s observed at radio frequencies. Data taken from this paper ; McLean et al. (2012*): B erger et d] 
<201(]h : lRoute & Wolszcza3 ( l2012r ): [Stelzer et alJ ( [20T^ : lGrosso et"^ d2007l) : ISchmitt & Liefk^ ( 120041) : IBurgasser et al.ii2003a) : iFleming et aUdigOal) 
and references therein. 
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Spectral type 



Fig. A.l. Ha luminosity - spectral type for all ultracool dwarfs observed at radio frequencies. Squares represent luminosities w hile 
triangles show upper limits. Data are taken from the Uterature dMcLean et al.1120121: iBurgasser et al.]|2003at feerger et al.ll2010l and 
references therein). 
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Fig.A.2. X-ray luminosity - spectral type for all ultracool dwarfs ob served at radio frequencies. Circles represent luminosities 
while triangles show upper limits. Data are taken from the literature (iMcLean et al. 2012t IStelzer et"ani2012t iBerger et"al]|2010l : 
iGrosso et al.1l2007t ISchmitt & Liefkell2004t iBurgasser et"ani2003at [Fleming et al.lll993L and references therein). 



